      function tlake(ix,iy,timeho)
      use alloc
      implicit real*8(a-h,o-z)

      pi180=4.*datan(1.d0)/180.

      tlake=tsmed(ix,iy)

      return
      end
